NASA TECHNICAL 


CNI 

CO 

CN| 

I 

X 


un 


MEM ORAN D U M 



>LOAN COPY: RETURI IJj 
AFVvL(DOGL) ~ 

KIRTLAisiD AFB, N. iVi. 




SERT II SPACECRAFT 
ATTITUDE RESPONSE 


by Bruce E. LeRoy, John D. Regetz, 
and Robert R. Lovell 

Lewis Research Center 
Cleveland, Ohio 44135 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • JULY 1971 


OBbTSTO 


TECH LIBRARY KAFB, NM 

iniiiiini 


1. Report No. 2. Government Accession No. 

NASA TM X-2324 

U JJ J JU IOU 

o, necipiem s uaraiog mo. 

4. Title and Subtitle 

SERT II SPACECRAFT ATTITUDE RESPONSE 

5. Report Date 

July 1971 

6. Performing Organization Code 

7. Author(s) 

Bruce E. LeRoy, John D. Regetz, and Robert R. Lovell 

8. Performing Organization Report No. 

E-6218 

10. Work Unit No. 
120-26 

9. Performing Organization Name and Address 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 

16. Abstract 

This report presents an analysis of the SERT II attitude response. Detailed analyses of magnetic, 
solar pressure, and ion thruster disturbance torques are presented. The attitude responses indi- 
cated by linearized equations of motion and a complete simulation of the satellite are compared to 


flight data. The discrepancies between flight and predicted data are attributed to errors in the 
estimation of satellite magnetic moments. The use of magnetic moments derived from the flight 
data resulted in a good correlation. Thrust vector orientation of the ion engines was a potential 
problem before the flight, so the capability of gimbaling the ion engines was provided. In flight, 
however, the misalinements were found to be small and no corrections of the thrust vector orien- 
tations were made. 


17. Key Words (Suggested by Author (s) ) 

Attitude control Magnetic moments 

Control moment gyros Ion thrusters 

Gravity-gradient 


18. Distribution Statement 

Unclassified - unlimited 


19. Security Classif. (of this report) 
Unclassified 


20. Security Classif. (of this page) 

Unclassified 


21. No. of Pages 
66 


22. Price* 

$3.00 


For sale by the National Technical Information Service. Springfield. Virginia 22151 



SERT II SPACECRAFT ATTITUDE RESPONSE 


by Bruce E. LeRoy, John D. Regetz, 
and Robert R. Lovell 

Lewis Research Center 
SUMMARY 

This report presents an analysis of the SERT II attitude response. Detailed anal- 
yses of magnetic, solar pressure, and ion thruster disturbance torques are presented. 
The attitude responses indicated by linearized equations of motion and a complete simu- 
lation of the satellite are compared to flight data. The discrepancies between flight and 
predicted data are attributed to errors in the estimation of satellite magnetic moments. 
The use of magnetic moments derived from the flight data resulted in a good correlation. 
Thrust vector orientation of the ion engines was a potential problem before the flight, 
so the capability of gimbaling the ion engines was provided. In flight, however, the mis- 
alinements were found to be small and no corrections of the thrust vector orientations 
were made. 


INTRODUCTION 

The SERT II spacecraft (Space Electric Rocket Test) was launched on February 3, 
1970, into a near-polar "sun-synchronous” orbit. The primary mission objectives 
were a 6 -month life test of a mercury -bombardment ion thruster system and the deter- 
mination of the thruster system operating characteristics in a space environment. 

Power for the experiments and housekeeping was provided by two foldout solar panels, 
which are rigidly fixed to the body. 

The orbiting vehicle consists of an Agena with the solar arrays mounted to the aft 
equipment rack and the spacecraft and spacecraft support unit mounted to the forward 
equipment rack. The inertial distribution of this configuration is such that it is gravity- 
gradient stabilized. The axis of minimum moment of inertia normally alines with the 
local vertical, and the axis of maximum moment of inertia alines with the orbit plane 
perpendicular. Two control moment gyroscopes (CMG's) are added to provide rate 


damping. After orbit injection, the spacecraft is initially oriented by the Agena cold 
gas attitude control system. 

This report presents an analysis of the SERT n attitude response. Preflight esti- 
mates of attitude perturbations and the resulting satellite response are compared to 
actual flight data. Based on the flight data, the satellite magnetic moments are calcu- 
lated and compared to measurements taken on the ground. The satellite response to ion 
engine misalinements is also analyzed and estimates of the misalinements obtained. 

The report derives in detail the solar pressure, magnetic, and ion thruster mis- 
alinement torques. Although gravity-gradient control moment gyro stabilization is not 
new (refs. 1 to 4), the report derives the satellite equations in a form different from 
those in the references. 


CONFIGURATION AND AXIS SYSTEMS 

The orbiting vehicle consists of an Agena with solar arrays mounted to the aft 
equipment rack and the spacecraft and spacecraft support unit (SSU) mounted to the for- 
ward equipment rack. As shown in figure 1, the ion thrusters are mounted to the space- 
craft and alined so that their thrust vectors pass nominally through the vehicle center 
of mass. Two control moment gyroscopes (CMG's) are mounted in the spacecraft sup- 
port unit. 

Three body-centered axis systems are used in this report. They are the reference 
axis system, the geometric axis system, and the body (or principal) axis system. The 
reference axis system is defined as follows: the axis is perpendicular to the earth- 
satellite radius vector, in the orbit plane, and positive in the direction of the satellite 
velocity vector; the axis is perpendicular to the orbit plane and opposite in sense 
from the orbit angular momentum vector; the axis is in the orbit plane and oppo- 
site in sense from the radius vector. Figure 2 shows the relation between this coor- 
dinate system and the earth. The geometric axis system is fixed in the satellite and is 
defined as follows: the Xq axis is perpendicular to the longitudinal axis of the Agena 
and parallel to the longest edge of the solar arrays and, when the satellite is properly 
oriented, is positive in the direction of the orbit velocity; the Yq axis is perpendicular 
to the longitudinal axis of the Agena and perpendicular to the plane of the solar arrays; 
the Zq axis is along the longitudinal axis of the Agena and is positive in the direction 
of the spacecraft support unit. Figure 2 also shows the satellite orientation when the 
two coordinate systems are coincident. The body (or principal) axis system is defined 
by the principal moments of inertia. The x axis (roll axis) is the axis of intermediate 
moment of inertia. The y axis (pitch axis) is the axis of maximum moment of inertia, 
and the z axis (yaw axis) is the axis of minimum moment of inertia. Figure 3 shows 
that the body and geometric coordinate systems are nearly coincident. Using the com- 
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puter program described in reference 5, the direction cosine matrix relating the body 
and geometric systems is found to be 
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■ 1. 000 
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0. 0024' 
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-0. 0024 

-0.0080 

1. 0000. 
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(All symbols are defined in appendix F. ) 

The principal (or body) coordinate system and the reference system are related by 
three Euler angles (ip, 9, ep). As shown in figure 4, the yaw angle ip is a rotation 
about Zpj, the rotation about y' is a pitch angle 9, and the rotation about x" is a 
roll angle <p. These three angles are used to define the satellite attitude with respect 
to the reference coordinate system. The two coordinate systems are related by 


"x" 


Cxp C 9 

Sip C 9 

-S 9 

"x" 

y 

= 

Cip S9 Sep - Sip Cep 

Cip Cep + S9 Sep Sip 

C9 Sep 

Y 



Cip S0 Cep + Sp Sep 

Sp S0 Cep - Ci p Sep 

C 9 Cep. 

JZ_ 


( 2 ) 


The orientation of CMG axes with respect to the body axes is shown in figure 5. 
Electromagnetic torquers are used to maintain the spin reference axes in the "V" con- 
figuration ±28° from the negative pitch axis. 

Table I gives the SERT II vehicle and orbit parameters, and table n gives the CMG 
parameters. 


SATELLITE DISTURBANCE TORQUES 

A complete analysis of solar pressure, magnetic, and ion thruster misalinement 
torques is given in appendix A. For a linearized analysis, we note here the estimated 
solar pressure and magnetic torques. 


Solar Pressure 


In appendix A, the solar pressure torque is calculated in two parts. Torque from 
solar pressure on the solar array is considered separately from the torque due to solar 
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pressure on the Agena- spacecraft. Under some simplifying assumptions (small Euler 
angles, coincidence of the geometric and principal axes), it was found that the torque 
from solar pressure on the Agena- spacecraft was an order of magnitude smaller than 
the torque due to the solar array. Since the torque from the Agena-spacecraft tends to 
cancel the torque from the solar array, a "worst case" analysis may be made by con- 
sidering only the torque from the solar array. In component form, the solar array 
torque is given below. 

The estimated maximum roll torque occurs when the apparent sun declination 5 is 
zero: 


T x = _3 ’ 40x10 4 N_m (3) 

The estimated maximum constant pitch torque occurs when 6 = 0 and the yaw angle 
is iO. 1 radian: 


T y = ±2.26xl0“ 5 N-m (4) 

The estimated sinusoidal pitch torque of maximum magnitude occurs when 6 = 30° (the 
maximum apparent declination): 


T y = 9. 82x10“ 5 sin wjt 


(5) 


Magnetics 


The estimated magnetic torques are found from the approximated earth field given 
in appendix A and the estimated satellite dipole given in appendix B. 

The approximate earth field in the body axis system given by equation (A31) is 


H = 


0. 1965 cos co*t 
0.03 


0. 3915 sin co*t 


oersteds 


( 6 ) 


From appendix B the worst- case satellite dipole moments are 

\4n 


P = 


-9.8x10' 
-8.2X10" 
.-9. 8xl0 4 J 


A-cm^ 
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Then from equation (A30) the estimated magnetic torque is 

[-3. 21X10 3 sin co*t + 2. 94xl0 2 1 


T B = 


4. 29xl(r cos co*t dyne-cm 

1. 611X10 3 cos wjt - 2. 94X10 2 


-3. 21X10" 4 sin co*t + 2. 94xl0~ 5 

T d = 4. 29X10' 4 cos w *t N-m 0 

D O 

.1. 61X10 -4 cos w*t - 2. 94x10" 3 

Equations (3), (4), (5), and (7) are used to estimate the satellite attitude response. 


SATELLITE ATTITUDE CONTROL TORQUES 

A complete analysis of gravity -gradient CMG control is given in appendix C. The 
linearized control torques in the body (principal) axis system are 

Linearized gravity-gradient torques: 

Tggx = -3K> 2 «yy - y? < 8a 

T eE , - - W 9 < 8b 


T =0 
ggz 


Linearized Euler cross-coupling torques: 


T ccx = -K) v " W* 


T =0 
ccy 


T ccz = "K>\y " X xx^ 
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Linearized CMG torques: 


T CMGX = cos 28° + 2a cos 28° + a>* cos 28°(aj + Og) J (10a) 

T CMGY = sin 28° [20 sin 28° - a>*(cjj - crg)cos 28° - 4e sin 28° cos co*tJ (10b) 
2 

T CMGZ = “ — cos^28°fa)*(aj + Og) + + 2a)*(^'J + 2h<p cos 28° - 2hu>*i// cos 28° 

C (10c) 

Equations (8) and (9) show that restoring torque exists about all three principal 
axes. Equation (10) shows that the CMG’s directly damp pitch and yaw rates and pro- 
vide additional restoring torque about the yaw axis. Equation (10a) is misleading as 
shown for it appears that the CMG's provide restoring torque about the roll axis. A 
transient restoring torque in roll exists due only to CMG lags. If it were possible to 
maintain a constant roll offset while holding pitch and yaw offsets to zero, each CMG 
gimbal angle would attain a steady-state position of -cp. Thus the second and third 
terms of the equation cancel. Although roll rates are not damped directly, the roll axis 
is the axis of intermediate moment of inertia, and roll energy will be cross coupled 
into the pitch and yaw axes. Thus, complete three-axis control is accomplished with 
gravity-gradient CMG torques. 


SPACECRAFT RES PONSE 
Linearized Equations of Motion 

Pitch . - Appendix C derives the linearized equations of rotational motion of the 
SERT II satellite. From equation (C28) we see that only pitch motion is affected by 
eccentricity. Using the parameters in table I, we find the pitch eccentricity response 
to be sinusoidal with a zero-to-peak magnitude of 

1 0 | = 5. 04 e for c = 2. 25 (N)(m)(sec) (Ha) 

1 0 | = 4. 16 e for c = 4. 5 (N)(m)(sec) (lib) 

Equations (C28) are also used to find the steady- state pitch response to constant and 
sinusoidal disturbances. Assuming e = 0, the steady- state step torque response is 
given by 
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( 12 ) 




where T y 
given by 


step 


= 56. 6 T 


- W 


is in newton-meters. The steady-state frequency response magnitude is 


9 SS | =197.8|T y | 

for c = 2. 25 (N)(m)(sec) 

(13a) 

Ksl = 172. 7|T y | 

for c = 4. 5 (N)(m)(sec) 

(13b) 


Roll . - Equations (C29) are used to find the steady- state roll response to step and 
sinusoidal disturbance. The step response is given by 


‘'’step = ^ 30 - 2 T x <»«> 

\y - W 

where T x is in newton-meters. The steady- state frequency response magnitude is 
given by 


Ks| t = 40. 3|T X | 

(15) 

«> S s| T =-i 40 - 3 lTj 

1 1 Z 

(16) 


where T„ and T_ are in newton-meters. 

X Z 

Equations (15) and (16) show that the roll response to sinusoidal roll torque is 90° 
out of phase with the roll response to sinusoidal yaw torques. 

Yaw . - Equations (C29) are also used to find the steady- state yaw response to step 
and sinusoidal disturbances. The step response is given by 


* 


^ K> 2 «yy - W - 2h < cos 28 ° 


= 96. 9 T 


(17) 
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where T is in newton-meters. The steady- state frequency response magnitude is 
given by 


* 


ss 


lb 

Y SS 


= 286. 3 T 
T„ 1 z 


= 207. IT 


j 40. 3 [ T x | 

for c = 2. 25 (N)(m)(sec) 
for c = 4. 5 (N)(m)(sec) 


( 18 ) 


(19) 


where T v and T are in newton-meters. As in the roll case the yaw response to T 
x z x 

is 90 out of phase with the yaw response to T . 


Results from Linearized Equations of Motion 

Equations (11) to (19) may be used to estimate the satellite response to various 
disturbances. 

Eccentricity . - Using equation (11), we find the maximum response magnitude for 
the worst-case eccentricity of 0. 005 is 

1 0 | = 0.0252 radian = 1.44° 

Solar pressure. - Equations (3) to (5) are used to estimate the solar pressure ef- 
fects. From equations (14) and (3), the maximum roll offset from solar pressure is 

cp m ov = -1. 027X10" 2 radian = -0. 59° 
r max 

From equations (12) and (4), the maximum constant pitch offset is 

9 = ±1. 279X10" 3 radian = ±0. 07° 

From equations (13) and (5), the maximum pitch oscillation is 

1 9 | =1. 94x10" 2 radian = 1. 11° 

Geometric axis misalinement . - Appendix C shows that the attitude is stable and 
will aline the principal axes with the reference axes. The horizon sensors are alined 
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with the geometric axes and will record a roll and pitch offset when the principal axes 
are alined properly. From equation (1), we see that these offsets are 

cp = -0. 008 radian = -0. 46° 

0 = 0. 0024 radian = 0. 14° 

if/ = +0. 06° 


Magnetics . - Equation (7) is used to estimate magnetic effects. Using equations 
(12) to (19), the satellite response to magnetic torques is 


«„l = 4 - 86< 


k s tepl=°- 05 ° 


- 0. 74 

7-ss T 


l^stepl = °- 16 ° 

M Tx -°- 740i 

Mt ,-*- 640 


Total effects . - Assuming all pitch sinusoidal disturbances are in phase, the maxi- 
mum pitch frequency response amplitude is 


ss 


max 


7. 41° 


The maximum constant pitch offset seen by the horizon sensors consists of geometric 
axes misalinements and solar pressure effects 



max 


0 . 21 ° 
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L 


The maximum roll oscillation magnitude is 


The maximum constant roll offset seen by the horizon sensors consists of geometric 
axes misalinements, solar pressure, and magnetic effects. 

|^step| max 

The maximum yaw oscillation magnitude is 


* 


ss 


max 


2. 74° 


The maximum constant yaw offset is 


* 


step 


max 


0 . 22 ° 


The response indicated above is the worst-case response found using the linearized 
equations of motion. The pitch frequency response magnitude appears to be the one item 
worthy of concern. If the magnitude is as given, only a small steady- state offset is 
allowed before the pitch angle exceeds the linear range of the horizon sensors (±10°). 
This is particularly important if the ion engines are initially misalined. However, the 
simulated response, described in the following section, showed that the pitch oscillation 
magnitude is approximately 4°, a tolerable value. 


Predicted Response 

Appendix D briefly describes the analog and digital simulations of the satellite dy- 
namics. If it is assumed that the satellite dipole is 


P B “ 


-9. 8X10 


-8. 2X10 


-9.8x10' 


A-cm' 


and that the CMG gain is h/c = 2. 0, the satellite steady-state attitude response was 
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simulated. Figure 6 shows the worst-case predicted response compared to actual flight 
data. 

The magnitudes of the simulated attitude and CMG angles are all larger than the 
flight data. The simulation graph is keyed on the pitch response. So with the simulated 
pitch response in phase with the flight data, all other simulated responses are out of 
phase with the flight data. 

The flight data shown in figures 6 and 7 are from one orbit during the first month 
of the mission. The data do not repeat exactly from one orbit to another, but the differ- 
ences are minor. The magnitude of the pitch oscillation is the parameter which varies 
the most. The variation is less than 1°, and this may be attributed to a changing ori- 
entation of the earth's magnetic dipole. The response repeats over a 24-hour period. 

The flight data are shown in steps since the parameters are sampled, and the values 
held for 1 -minute intervals. In addition, the telemetry system quantizes the data so that 
pitch readings from the horizon sensors appear at 0. 5° steps with the associated reso- 
lution of ±0. 25°. The roll steps are 0. 4° and the roll resolution is =t0. 2°. Since the 
CMG gimbal angle variation is large (~±30°) and the allowed number of discrete steps 
is fixed (~60), the gimbal angle step size is approximately 1 . 20 ° and the resolution 
is 0 . 6 °. 


RESULTS FROM FLIGHT DATA 
Magnetic Moment Estimation 


Estimation of the satellite magnetic moment could be easily accomplished if flight 
data on all three Euler angles were available. Since direct measurement of the yaw 
angle is not possible, CMG information is used. From equations (C29), the steady- 
state frequency response magnitude of the CMG gimbal angle sum is given by 


a" = 232.6 T for c = 2. 25 (N)(m)(sec) 
z 


116. 6 T z f or c _ 4 # 5 (N)(m)(sec) 


( 20 ) 


where <x" = a 1 + <j 0 and T is in newton-meters. 

12 z 

If we assume that the earth's magnetic field in the body axis system is given by 
equation ( 6 ), equation (A30) may be approximated by 
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dyne -cm 


( 21 ) 


p y (0. 03915 sin o>*t) 



p (0. 01965 cos «*t) - p„( 0. 03915 sin a>*t) 

Z O A U 


-p y (0. 01965 cos u>*t) 


— o 

where p is in ampere-cm and the terms containing Hy are neglected. From equa- 
tions (15) and (16) the roll frequency response magnitude is given by 


I <P 


ss 



( 22 ) 


where T and T are in newton-meters. 

X z 

Although the expected eccentricity was 0. 005, the achieved eccentricity was 0. 0003, 
hence, the pitch oscillation magnitude is assumed to be a result of magnetics and solar 
pressure. Ground test of CMG's showed a high probability of CMG operation near a 
gain of 2. Thus, assuming the CMG gain is 2.0, equations (13), (20), and (22) maybe 
solved for the magnitude of the disturbance torques: 


| Ty | = 5. 056xl0~ 3 \ d | N-m 

(23a) 

|T Z | = 4. 3xl0 -3 |a" | N-m 

(23b) 

|T X | 2 = 6. 16xl0" 3 |<p j 2 - 1 T z | 2 N-m 
Using the flight data and solving equations (23) gives 

(23c) 

| T v | =2. 38X10" 4 N-m 

(24a) 

|Tj = 9.53X10" 5 N-m 

(24b) 

|T X | = 1. 97xl0 -4 N-m 

(24c) 


From equation (21), 


|T X | = 0. 03915 |p y | 
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and 


|Tj = 0. 01965 |p y | 

The average of the two estimates of p y is 

|p y | = 4. 94X10 4 A-cm^ 

The flight data shown in figures 6 and 7 were recorded when the sunline deviation 

from orbit normal was large. Thus, the pitch response is subject to both solar and 

magnetic torques. Since the sunline deviation was approximately 30°, the pitch torque 

-5 

due to solar pressure was found to be 8. 85x10 newton-meter. Assuming the solar 
and magnetic torques are in phase, the pitch magnetic torque is 

| T y | = 1. 495X10" 4 N-m 


Now assuming p - p = p, 

X Z 


|T | = 0.0438 p 

and the initial estimates of |p x | and |p z | are 

|p x | = |p z | = 3. 41X10 4 A-cra 2 


The polarity of the magnetic moments may be determined by observing the roll and 
pitch angles when the satellite is near the north magnetic pole and the magnetic equator. 
The resulting estimated total dipole is 


P B ~ 


3. 41X10 


-4. 94x10 


-3.41X10 


A-cra' 


(25) 


Figure 7 shows the simulated response using the estimated dipole. The figure 
shows good agreement of roll and pitch angles between the simulated and flight data. 

The simulated cr^ response is of the proper magnitude and phase but is offset from the 
flight data. The offset may be due to an incorrect torquer torque level. The torquers 
are not regulated but are subject to temperature and current variations. The simulated 
o ' 2 response shows the CMG input rates are sufficiently small that stiction prevents 
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motion. The simulated Og response is offset from the flight data, but its offset is not 
as large as the offset. 

The dipole in equation (25) is a good approximation to the true dipole. Varying the 
roll and yaw magnetic moments will change the simulation from that shown in figure 7. 
However, these changes are of such a nature that the approximation to flight data may 
not be considered improved. 

Using the spacecraft/SSU dipole given in appendix B, the resulting Agena/solar ar- 
ray dipole is 



_ 3. 41X10 4 + 2. 27X10 4 


5. 68X10 4 

p = 

-4. 94X10 4 + 1. 11X10 4 

= 

-3.83X10 4 


.-3. 41X10 4 + 2. 33X10 4 . 


-1.08X10 4 . 


ion Engine Response 

The simulations described in appendix D were used to determine acceptable ion en- 
gine thrust misalinements. It was found that if | a\ and/or |/ 3 | were approximately 4° 
with the engine operating at full thrust, the satellite's step response stability limit was 
exceeded and the satellite would tumble. To avoid this problem, provisions were made 
to allow the engines to operate at 30, 80, and 100 percent thrust levels. In addition, 
each engine was mounted on a gimbal system capable of correcting a thrust error as 
large as ±10° in the a or /3 axes. Thus corrections were to be made to each thrust 
level to avoid exceeding the stability limit. 

It was found, however, that the thrust vectors were reasonably well alined and no 
gimbaling of the engines was performed. If the ion thrust misalinement is small (a ^ 10° 
and /3 < 10°), the roll and pitch misalinement torques given in equation (A38) may be 
approximated by 


T x = -l z F sin /3 


(26a) 


T = -F 

y 


V 


2 2 

l v + l _ sin a cos /3 

X Z 


(26b) 


If the Euler angles 9 and cp are also small, equations (12), (13), and (26) maybe used 
to find the thrust misalinements. 

Figure 8 shows the steady- state roll and pitch horizon sensor readings before and 
after the thrusters were started. Before the thrusters were started, the steady- state 
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offsets were 


9 = +0. 20° 

(27a) 

cp = -0. 65° 

(27b) 


The offsets were found by averaging the horizon sensor data over the orbit. After the 
fore thruster was started and operating at 100 percent thrust, the horizon sensor read- 
ings showed steady- state offsets of 


9 = -1. 96° 

(28a) 

o 

m 

00 

H 

II 

(28b) 


The offsets given in equations (27) are assumed to be the satellite response to principal 
axis misalinements, solar pressure, and magnetics. Thus the roll and pitch response 
to the thruster is 


9 = -2. 16° (29a) 

cp = - 1.20° (29b) 

Using equations (29), (26), and (14) and the parameters in table II gives 

0 = 0. 51° (30) 

Now using equations (12), (26), (29), and (30) and solving for a gives 

a = 0. 48° (31) 

After the aft thruster was started and operating at 100 percent thrust, the horizon sensor 
readings showed steady- state offsets of 

9 = +1. 20° (32a) 

cp = -0. 35° (32b) 
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Following a similar procedure gives the misalinements as 


13 = -0. 13° 

(33a) 

a = -0. 22° 

(33b) 


CONCLUDING REMARKS 

This report has presented an analysis of the SERT II attitude response. The mag- 
netic, solar pressure, and ion thruster torques were examined in detail. The analysis 
and simulation predicted an attitude response tolerable to mission objectives and did not 
reveal any instabilities under normal operation. The analysis, however, did show a 
potential instability due to thruster misalinement. The potential instability was elim- 
inated by providing gimbaling for the ion thrusters. 

Flight data compared to the predicted response showed that the analysis was correct 
within the tolerances on such parameters as magnetic moments, inertias, CMG gains, 
etc. Flight data were used to estimate the magnetic moments. The simulated response 
using the estimated moments showed good agreement with the flight data. Further im- 
provement of flight and simulation data agreement is unlikely unless the earth’s mag- 
netic field model is improved and the CMG parameters are more accurately known. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 6, 1971, 

120-26. 
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APPENDIX A 



SATELLITE ENVIRONMENT AND DISTURBANCES 

To estimate or model the satellite attitude response, it is necessary to analyze the 
spacecraft environment and possible attitude disturbances. It was found that solar pres- 
sure torques, magnetic torques, torques from principal axes - geometric axes misaline- 
ments, and ion thruster misalinement torques would produce attitude response detect- 
able by the onboard horizon sensors. Other sources of attitude errors were found to be 
negligible. 


Solar Pressure 

Derivation of torque equations . - For purposes of solar pressure analysis, the ve- 
hicle is considered to be a cylinder and a plate. The aft equipment rack and spacecraft 
geometric irregularities are ignored. References 6 and 7 give the differential force on 
the differential area to be 

dF = |(1 + R)(s • n) 2 n + (1 - R)(s • n)£n x (s x n)J^(- P g )dA (Al) 

The differential torque is given by 


dT = r B X dF (A2) 

The apparent sun motion as viewed from the spacecraft is a circle if the body angles 
are zero. We define an Xg, Yg, Zg system coincident with the reference axis system 
at some point in the orbit such that a simple rotation about the Xg axis places the new 
Y axis (Yg) along s. This is illustrated in figure 9. The Xg, Yg, Zg and Xg, Yg, 

Zg systems are related by 


'X' 

T 

1 

o 

o 

i 


"X' 

Y 

= 

0 cos 6 -sin 6 


Y 

_Z. 

B 

.0 sin 6 cos 6. 


_Z_ 


(A3) 
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The reference system and the Xg, Yg, Zg system are related by 


*x‘ 


'cos(u) 0 t + y) 

0 

sin(o> o t + y)' 

'x' 

t 

Y 

= 

0 

1 

0 

Y 

(A4) 

.Z. 

R 

.-sin(o> o t + y) 

0 

cos(co Q t + rl 

.Z. 

B 


where y is used as a phase angle to relate the systems in time. 
The unit vector s in the S system is given by 

'0 

s = 1 
-0J 

Thus, in the reference system, the vector s becomes 


*"sin 6 sin(w Q t + y) 


S R = 


cos 5 

sin 6 cos(w t + y) 


By using equation (2) and assuming that the body rotation angles are small, a vector in 
the reference system may be found in the body system by 


*1 i p 

1 

Le -<p 

Thus, in the body system, the vector s becomes 


-9 

V 

1 J 


rx 

Y 

ZJ 


R 


S B " 


S6 S(a) Q t + y) + ^ C6 - 9 S6 C(u Q t + y) 
S6 S(oj Q t + y) + C6 + cp S6 C(u> 0 t + y) 
9 S6 S(cu Q t + y) - <p C6 + SS C(cu 0 t + y) 


(A5) 
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For the solar panels, the vector n in the body system is 


n B “ 


- 0. 0010 
1. 0000 
-0. 0080. 


Thus, Sg X Sg = -0. 001 S v + S v - 0. 008 S 7 and 


x y 


n B x (s B x n B ) = 


S v (l + 6. 4xl0“ 5 ) + lxlO" 3 S v - 8xl0" 6 S, 
x y z 

1x10“ 3 S v + 6. 5X10 -5 S„ + 8xl0 -3 S„ 
x y z 

-8X10' 6 S v + 8X10" 3 S, r + S(1 + 10" 6 ) 
x y z 

where S , S , S are given by equation (A5). Substituting in equation (Al) and integrat- 
x y z 

ing over the area gives 


f B = - p s A(S y ” °- 001 S x " °- 008 S 7 )^ 


(1 + R)(S y - 0. 001 S x - 0. 008 S z )(-0. 001) 

+ (1 - R) [s x (l + 6. 4xl0 -5 ) + lxlO -3 Sy - 8x10”® sj 

(1 + R)(S y - 0. 001 S x - 0. 008 S z ) 

+ (1 - R)(6. 5x10”® S + lxlO' 3 S + 8X10' 3 S ) 
y x z 7 

(1 + R)(S y - 0. 001 S x - 0. 008 S z )(-0. 008) 

+ (1 - R)[s z (l + 10' 6 ) + 8xl0” 3 S y - 8x10'® S J 


(A6) 


The solar panels are offset in the +Y geometric axis direction. Thus, the radius 
vector to the panels’ center of pressure in the geometric system is 


r G " 


0 


GY 


l GZ 
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Using equation (1) the radius vector in the body system is 



-0. 001 Tqy + 0* 0024 r GZ 


r x 

r B " 

r Gy + 0. 008 r GZ 

= 

r y 


-0. 008 r GY + r GZ 


_ r z_ 


and the torque on the body due to solar pressure on the solar panels is 

F B = *< r y F z ' r z F y> + y (r z F x " r x F z ) + £(r x F y “ r y F x> < A7 > 

To determine the solar pressure forces on the Agena and the spacecraft, we intro- 

a . s. 

duce an additional coordinate system. The i, j, k system shown in figure 10 is a right- 
handed coordinate system with the origin at the center of mass and with k along the 
axis of the cylinder, i directed away from s such that the 1-k plane contains s, and 
f completes the set such that j • s = 0. The angles a' and >3 T are defined in fig- 
ure 10. 

A A X 

In the i, j, k system 


n = 


-sin a' 
-cos a ’ 
0 


and 


-cos /3' 
0 

-sin /3' 


The differential area is given by dA = r* d a' dh\ Thus, equation (Al) becomes 
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dF = -P s < 

r 

(1 + R)(cos (3 T sin 

a') 2 

-sin a' 
-cos a' 




0 


+ (1 - R)(cos 0' sin ex') 


2 

-cos a ’ cos 0' 


lr’ do’ dh’ 


| sin a ' cos a?’ cos 0* 
-sin 0' 

The moment arm to the differential area in the i, j, k system is 

-r' sin a ' 
rg = | -r' cos a ' 


Thus, the differential torque from the differential area is 
dT = i(-r' cos a ' dF^ - r^ dFj) + ^(r^ dF^ + r' sin a-' dF^) 

+ k(-r' sin a' dFj + r' cos a ' dF.) 

where the components of dF are given in equation (A8). 

Integrating over the area exposed to sunlight gives the torque in the i, j, k 

T = i(0) + fjp g r’ cos 0’h’ jr fc cos 0^2 + ^ rJ + 1 r’(l - R)sin 0'jj + k(0) 


Now, a vector in the i, j, k system may be found in the geometric system 




(A8) 


(A9) 

system: 

(A10) 
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i z X 


-sin I tan 


-1 s gx n 
3 gy/ 


-cos (tan 


■1 S GX^ 


S, 


-cos (tan -1 — ^ 




3 gy> 


sin [tan 

\ S 


GY/ 


1 s gx n 


GY/ 


-lJLkJ 


(AH) 


where Sq X and S^y are the x and Y components of s in the geometric system. 
Using equation (1) we find 


S GX = s x + °* 001 s y " °- 0024 s z 
S„ v = S„ - 0. 001 S Y - 0. 008 

Lr jl y X Z 

where S x , Sy, S z are the components of s in the body system and given in equa- 
tion (A5). 

Using equations (A10) and (All) we find the torque in the geometric system: 



cos 


/3 r h' jr k cos 0^2 + | rJ + 1 r'(l - R)sin 


cos /3’h'Jr k cos /S’ ^2+^ R ^ +~ r ' U ~ R)sin /3'Jj^ + Zq(0) 

(A12) 


where 


sin /3’ = Z n • Sp. = S + 0. 0024 S v + 0. 008 S„ 
kj z x y 


and 


cos /3' = yl - sin^/3’ 

Using equations (1) and (A12) we find the torque in the body system to be 
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t B = x(T GX - 0. 001 T gy ) + y(T GY + 0. 001 T gx ) + z(-0. 0024 T gx - 0. 008 T QY ) 

(A13) 


Equations (A7) and (A13) give the total solar pressure torque on the satellite. 

Approximate torque equations. - Estimation of the solar pressure torque is not 
easily accomplished with the equations given. To estimate the torque from solar pres- 
sure on the solar array we first assume that n B is given by 



then 


s 


B 



and 


n B X X ”B^ “ 


x 

0 


- S z. 


The components S v and S„ are given in equation (A5). Using these approximations, 

X z 

equation (A6) become 


F B = - p s AS y 


(1 - R)S ? 

(1 + R)Sj 

(1 - R)S„ 


(A14) 


Since we limit the sun declination to 30° or less, |F | will always be larger than 
I F I or I F_ I . We may not, however, eliminate F„ or F from the approximation. 
Since the panels' reflectivity is small (R ~ 0. 2), when 6 ~ 30 u the force components are 
all the same order of magnitude. 
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For SERT H, the moment arm to the panels' center of pressure in the geometric 
system is 


r G = 


’ 0 


‘ 0 


• 0 * 

r GY 

= 

0. 20 

meters or 

0.67 

_ r GZ_ 


.-3. 35. 


.-11.0. 


ft 


In the body system, the moment arm becomes 


r B = 


-0. 008 
0. 18 
.-3.35J 

-0. 027 
0. 58 
-11.0 J 


or 


or 


r o 

0 

L-3. 35 

‘ 0 

0 

.-11 J 


meters 


ft 


Using the approximation to r^, equation (A7) becomes 


T B = X (" r Z F y) + y( r z F x ) + z(0) 


(A15) 


Noting that the Euler angles are small, we may approximate the force components 


as 


F x ~ -P S A(1 - R) jc6 S6 S(co o t + y) + # C 2 6] 
Fy ~ -PgA(l + R)C 2 6 

where we have retained the most significant terms. 

Letting 


(A16a) 

(A16b) 
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A = 17. 65 m 2 (190 ft 2 ) 


P c = 4. 788X10" 6 N/m 2 (1x10“ 7 lb/ft 2 ) 
s 

\p = i0. 1 
R = 0. 2 

we find that the maximum constant roll torque occurs when 6 = 0 and is 

T x = -3. 40X10" 4 N-m = -2. 51xl0“ 4 ft-lb 

The maximum constant pitch torque also occurs when 6 = 0 and xp = i0. 1 and is 

Ty = ±2. 26X10" 5 N-m = ±1. 67x10" 5 ft-lb 
The sinusoidal pitch torque of maximum magnitude occurs when 6 = 30° and is 

Ty = 9. 82x10"^ sin(co*t + y) N-m 

= 7. 24x10" 5 sin(a)*t + y) ft-lb 
The roll torque for 6 = 30° becomes 

T x = -2. 55xl0" 4 N-m 
= -1. 88X10" 4 ft-lb 

These expressions provide an estimate of the torque due to solar pressure on the solar 
panels. 

To estimate the solar pressure effects on the Agena and spacecraft, we first assume 
that the Euler angles may be neglected. Then the following relations apply: 
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S GX S x 

— r^i — - 

S GY S y 


= tan 6 sin(u>*t + y) 


-cos 




-cos 5 

V 1 - c 2 (w*t + y)s 2 6 


sin 




S5 S(co*t + y) 

^1 - S 2 6 C 2 (oj*t +y) 


sin j3' = S6 C(w*t + y) 

cos 0' = ifl - S 2 6 C 2 (w*t +y) 
We also assume that equation (A13) may be approximated by 



o 


' T x' 

t b = 

t gy 

1! 

T y 


J 

O 

1 


1 

o 

1 


Now letting 


r T = 0. 76 m (2. 5 ft) 

r k = -0. 36 m (-1. 17 ft) 

P o = 4. 788X10" 6 N/m 2 (lxlO -6 lb/ft 2 ) 
s 


R = 0. 8 


2r , h» = 9. 01 m 2 (97 ft 2 ) 

and using these relations, we find T and T from equation (A12). For 6 = 0, the 

x y 

torque due to solar pressure on the Agena-spacecraft is 
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T x ~ 1. 95X10" 5 N-m = 1. 44xl0” 5 ft-lb 

Ty ~ 0 

For 6 = 30°, the roll and pitch torques from solar pressure on the Agena-spacecraft can 
be approximated by 


T x = 1. 69X10' 5 N-m - 2. 24xl0' 6 cos(w*t + y) N-m 
= 1. 25X10" 5 ft-lb - 1. 65X10" 6 cos(ct>*t + y) ft-lb 

T y = -9. 76X10" 6 sin(w*t + y) N-m 
= -7. 2x10”® sin(w*t + y) ft-lb 

It is evident from the previous calculations that the torques from solar pressure on 
the Agena-spacecraft are an order of magnitude less than the torques from solar pres- 
sure on the solar array. In addition, for 6=0, the Agena-spacecraft roll torque is 
opposite in sign from the solar array roll torque. For 6 = 30°, the constant components 
of the roll torques are of opposite sign, while the pitch torques are 180° out of phase. 

The approximations indicate that solar pressure torques on the Agena-spacecraft tend to 
lessen the torques from the solar array. Thus, for crude analysis purposes, neglecting 
the Agena-spacecraft solar pressure torques would provide a worst case. 


Magnetics 

The instantaneous magnetic torque on the satellite is given by 


T = p x H 


(A17) 


where p is the satellite dipole and H is the earth field vector. If p is given in 
2 

amperes-cm and H is oersteds, the torque is given in dyne-cm by 

T = — (p x H) (A18) 

10 
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To find the torque in the body (principal) axis system, several coordinate definitions 
and transformations are needed. We define an earth-centered inertial coordinate system 
Xp Yp Zj with Xj pointing toward the vernal equinox, Zj along the earth- spin axis, 
and Yj completing the right-handed set. We also define an earth- centered "magnetic” 
coordinate system X^, Y M , Z M with Xj^ pointing toward the 0 ° longitude point on 
the equator, Z M positive northward and rotated 18° from the earth-spin axis about X^, 
and Y M completing the right-handed set. The two coordinate systems are related by 
the transformation: 


x" 


Cft E 

-sn E ci8° 

SS2 e S18°" 

"x" 


Y 

= 


Cn E C18° 

-Cft E S18° 

Y 

(A19) 

.Z_ 

I 

. 0 

S18° 

C18° 

.Z_ 

M 


where n E is the inertial longitude of the Greenwich meridian. 

Since the earth's field is an assumed dipole with axis Z^, we can express the field 
at the location of the satellite in terms of two components. The axial component is that 
component parallel to the dipole axis. The normal component is perpendicular to the 
axial component and directed from the satellite toward the earth's dipole axis. Figure 11 
shows the unit axial and normal vectors. When the unit vectors are defined in this man- 
ner, reference 8 gives the magnitude of the axial and normal components in oersteds as 


H = 0 . 308 1 - ~ 3 c ° s H 

aX q 

/ r c\ 


iR 


E/ 


H 


norm 


0 . 461 gH Lg S 



where 77 is the angle between the dipole axis and the radius vector to the satellite. The 
magnitude of H is 



+ H 

norm 
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We now define an earth-centered orbit system Xq, Yq, Zq such that Xq is along 
the line of nodes and positive toward the ascending node, Zq is positive along the orbit 
angular momentum vector, and Yq completes the right-handed set. The orbit and in- 
ertial systems are related by 


X 


CS2 Sfi 0 


X 

Y 

= 

-SS2 Ci cn Ci Si 


Y 

Z_ 

o 

. SS2 Si -Cfi Si Ci 


Z. 


(A21) 


Where fi is the inertial longitude of the orbit ascending node and i is the inclina- 
tion (i = 99. 1 for SERT n). The absence of a phase angle does not restrict the analysis. 
The inertial, magnetic, and orbit coordinate systems are shown in figure 12. The refer- 
ence axis system (defined previously) and the orbit system are related by 



(A22) 


The unit vector along the magnetic dipole axis is given in the magnetic coordinate 
system as 


H ax = 


Using equation (A19) H___ in the inertial system is found to be 

ax 


H 


ax 


SS2 Z S18° 
-CO, S18° 

Z i 

C18° 


Using equation (A21), IT 

ax 


in the orbit system is 


29 



H. 


ax 


S18° S(fi E - fi) 

-S18 0 Ci C(n E - to) + Si C18° 
S18° Si C(fi E - fi) + Ci C 18° 


Using equation (A23), H in the reference system is 

dX 


H ax = 


-Sw c t S18° S(n E - 0) + Cw 0 t(-S18° Ci C(fi E 
-S18 0 Si C(n E - O) - Ci C18° 
-Cw 0 t S18° S(n E - 0) - Sw 0 t(-S18° Ci C (fl E 


Q) + Si C18°) 
fi) + Si C18°) 


or 


H, 


ax 


H, 


Hr 


H, 


Now in the reference axis system, the unit radius vector from earth center to the 
spacecraft is given by 


0 


r 


R ~ 


0 


-1 


Thus, cos rj can be found by 

cos n = H ax • r R = -H 3 = Cu> 0 t S18° S(fi E - Si) + Sw 0 t(-S18° Ci C(« E - Si) + Si C18°) 

(A24) 

The unit normal vector H norm , as shown in figure 11, is perpendicular to and 

directed from the satellite to the dipole axis. Thus, H norm sin rj is given by 
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H, 


norm 


sin V = x (H« x r R ) 


ax R' 


(A25) 


and 


H 


sin r\ = 
norm 1 


~ H 1 H 3 

- H 2 H 3 
H? + H \ 


(A26) 


where Hp H 2 , Hg are given in equation (A23). 

Now, using equations (A20), (A23), (A24), and (A26), we find the earth's magnetic 
field in the reference system: 



0. 308 ^1 - SH^Hj + 0. 922(-Hg)(-H 1 Hg) 

- / R p \ 3 


Hr = ( n 

0. 308 (l - 3HgJH 2 + 0. 922(-Hg)(-H 2 Hg) 


0. 308^1 - 3Hg)Hg + 0. 922(-H 3 )(Hj + H 2 ) 


(A27) 


or 


H - I — 

h r L 

1 C> 


|H 2 (0. 308 - 0.002 HgJ 
Hg(0. 614 - 0.002 Hg) 

and since Hg max = 1 we may approximate H R by 


H 


i(°- 


308 - 0. 002 H 


H R = 


0. 308 Hj 1 
0. 308 Ho 


-0.614 H, 


(A28) 


For the SERT II expected orbit parameters given in table I, equation (A28) becomes 
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-0. 0614 Su> 0 t S(S2 e - a) + Ccv Q t (0. 00972 C(fi E - fi)+ 0. 1867)' 
-0. 0607 C(S2 e - ft) + 0. 0299 

0. 122 Cco Q t S(ft E - ft) + Soj Q t (0.0194 C(ft E - ft) + 0. 372) 


(A29) 


Now, if the Euler angles are small, Hg ~ Hg. If the satellite dipole is expressed 
as 





from equation (A18) the magnetic torque on the body is 


T = — 

1 -D 


B 


10 


Vz - Pz H y 

pA - PA 

^x H y ” ^y^x 


(A30) 


where H v , H , H_ are given by equation (A29). 
x y Zi 

For analysis purposes, equation (A29) may be approximated by 


0. 1965 cos o>*t 



0.03 


oersteds 


0. 3915 sin o>*t J 


(A31) 


where we have approximated the instantaneous orbit rate co Q by the average orbit rate 
u>*, let ft E - ft = 90°, and ignored the phase angles. 


Ion Thruster Misalinement 

Figure 13 shows a schematic of the geometric and thruster axis systems. The 
X T A’ y T A' Z T A s y stem is the reference system for the aft ion thruster. Similarly, 
the x T E , y T E , z T E system is the reference system for the fore thruster. In both 
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thruster reference systems, the x,j, axis is directed through the center of mass. The 
desired thrust direction is along x,p. 

The thrust misalinement is found in terms of the o and /3 angular errors defined 
by figure 13. The x”, y M , z" system is related to the x T , y T , z T system by: 


V 


'CaCjS 

-Cat S/3 

So* 

"x" 

It 

y 

= 

S/3 

C/3 

0 

y 

(A3 2) 

. z_ 

T 

.-So- C/3 

So S/3 

Co. 

_z. 



It is assumed that the true ion engine thrust is along the x M axis, and thus, 


F" = 


(A33) 


Using equation (A32) we find F in the thruster reference system: 


F T ' 


F Co- C/3 

F. S/3 

L-F S o C/3. 


(A3 4) 


Now the thruster reference system is related to the geometric system by 


x 


~ l * n 

l z 

Txl 




+ i z 


y 

= 

0 1 

0 

Y 

7 



~ l x 

7 

m ^ 

G 





(A3 5) 


where and l _ are defined in figure 13 and their values given in table n. Using 

A Z 

equation (A35), the thrust in the geometric system is 
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f g = 


-i x FCo C/3 - i z F So C/3 


V 


l l+ l l 


F S/3 


-L z F Co- Cj3 + l x F So- C/3 




,2 , 2 

*x + Z z 


(A36) 


The radius vector to the point of thrust application is given in the geometric system as 


r G " 


Thus the torque due to thrust misalinement in the geometric system is 


T G = 


F sin /3 


' ^ sin a cos & 


l x F sin /S 


(A37) 


Thus, using equation (1) we find the principal axis torques due to thrust misalinement 


are 


T B = 


•Z F sin /3 + 0. 001 F\ l ^ + Z ? sin a cos /3 + 0. 0024 FZ sin /3 

Z " X Z X 


0. 001 FZ sin /3 - fVz^ + Z^ sin o» cos /3 + 0. 008 Z F sin /3 

Z " X Z X 

[o. 0024 Z 2 F sin ,9- 0. 008 F ^ Z | + Z | sin a- cos 0 + Z x F sin 0 


(A38) 
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APPENDIX B 




MAGNETIC DIPOLE MOMENT MEASUREMENTS ON THE SERT II 
SPACECRAFT AND SSU 

To permit the determination of the magnetically caused disturbance torques acting 
on SERT II, it was necessary to obtain an estimate of the net magnetic dipole moment of 
the satellite in'its orbital configuration. This configuration can be separated into two 
parts for purposes of discussing the magnetic measurements which were made. The 
first part comprised the Agena stage and the solar array, for which estimates of the 
residual dipole moments were obtained from the Agena contractor. The second part con- 
sists of the spacecraft (S/C) and spacecraft support unit (SSU), which were either fabri- 
cated or integrated at Lewis and were, therefore, available for magnetic testing. The 
remainder of this discussion is concerned with measurements made on this latter part 
of the orbital configuration. 

Because the ion thrusters each contain several permanent magnets, it was expected 
that they would represent the largest single contributor to the net magnetic dipole mo- 
ment. An experimental setup was, therefore, devised to measure this dipole (see 
sketch a). The setup consisted of a large-diameter coil of copper wire of essentially 



zero length. The thruster was suspended by a torsion wire so that the center of the ef- 
fective magnetic dipole was as close to the center of the coil as possible. The entire 
setup was alined with the coil plane vertical and perpendicular to magnetic north. When 
a current of the proper polarity was passed through the coil, the period of oscillation of 
the suspended thruster was given by 


r 



I 

+ 10" 8 p T (H c + H e ) 


(Bl) 
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I III I 


where 

H e horizontal magnetic field intensity component of the earth, oersteds 

H c magnetic field intensity due to current at center of coil, oersteds 

o 

I moment of inertia of suspended assembly, kg-m 
K w torsional spring constant of suspension wire, N-m/rad 

2 

P T magnitude of magnetic dipole moment of thruster, A-cm 
With no current in the coil, the period is given by 



Combining equations (Bl) and (B2) gives 


Pq ’ — 


The value of was obtained by suspending a known moment of inertia from the 

wire and observing the period. This method was used to obtain the value of the dipole 

moment of a flight prototype thruster by suspending it in turn about each of the three 

principal axes. The dipole was found to be very nearly parallel to the thrust axis of the 

4 2 2 

thruster and to have a magnitude of about 23x10 ampere-cm . (One A-cm normal to a 

-8 

magnetic field intensity of 1 Oe has 10 N-m of torque exerted on it. ) 

Measurements on Overall Spacecraft-SSU 

Because of its large size and weight, methods which required dynamic motion of the 
S/C-SSU combination, such as the torsion wire method previously described, were felt 
to be impractical. Instead, it was decided that field measurement techniques would be 
used to obtain a value for the overall magnetic dipole moment of the combination. 

The magnetic field intensity along the axis of a magnetic dipole is given by the fol- 
lowing expression: 
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2pr 



where p is the magnitude of the magnetic dipole moment, r is the axial distance from 
the center of the dipole to the point of measurement, and l is the length of the dipole. 
If Z is assumed to be small compared to the measurement distance, a very simple ex- 
pression results for the magnetic intensity: 


H = 2£ 
r 3 

Thus, if the spacecraft is assumed to contain a single dipole of negligible length, an es- 
timate of the magnitude of that dipole can be obtained by measuring the magnetic inten- 
sity at a specified distance. This method was used to determine the dipole magnitude for 
the SERT n S/C-SSU combination. 

In the actual case, the S/C-SSU contains not one but several magnetic dipoles, each 
located in a different place within the combination. If the measurement distance is made 
sufficiently large, however, the dipole magnitude determined will be that of a single 
equivalent dipole. Since the torques produced by the interaction of the earth’s magnetic 
field with each of the actual dipoles are vectorially additive, determination of the mag- 
nitude of the single equivalent dipole is sufficient to establish the values of magnetic 
torque acting upon the satellite. 

The S/C-SSU combination was mounted in a fixture which allowed it to be both 
rotated about all three principal axes and translated horizontally. A sensitive magne- 
tometer was mounted on a tripod with the probe sensitive axis horizontal. The measure- 
ment procedure was as follows: With an S/C-SSU principal axis in line with the mag- 
netometer probe and as far from it as possible, the magnetometer was set to its most 
sensitive scale and any ambient field readings were biased out. The S/C-SSU combin- 
ation was then moved toward the probe until a specified field reading was obtained. The 
distance from the center of the combination to the center of the probe was measured. 

The combination was then moved away as far as possible, rotated through 180° to orient 
the negative of the same principal axis toward the probe, and the preceding steps re- 
peated. This procedure was repeated for the other two principal axes. Both positive and 
negative readings were taken on each principal axis to account for any offset of the cen- 
ter of the equivalent dipole from the center of the S/C-SSU combination. The positive 
and negative distance readings were averaged. 
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Results 


The torsional pendulum setup was used to determine the magnetic dipole moment of 

a flight prototype thruster, since these were expected to be the largest single contributor 

to the net dipole moment. The values obtained for this thruster are shown in table III. 

Also measured were several of the individual bar magnets which are used in the thrust - 

4 

ers themselves. The values obtained for these magnets averaged to about 3. 4x10 
o 

ampere-cm . As can be seen from the table, the main thruster dipole moment is es- 
sentially parallel to the thrust axis and has a magnitude of approximately 2. 3x10^ 
o 

ampere-cm . With the two thrusters mounted on the spacecraft in their nominal orien- 
tations, the thruster magnets have their N-poles alined toward the thruster exit on the 
orbit-raising thruster, while the N-poles are alined away from the thruster exit for the 

orbiting -lowering thruster. With this configuration, the net yaw-axis component is 

4 2 

zero, while the net roll-axis component is about 7. 6x10 ampere-cm . To reduce the 

A 

residual roll dipole to near zero, two thruster magnets having a total dipole of 6. 8x10 
ampere-cm , were attached to the S/C parallel to the roll axis. 

As a check on these results, the field measurement method was used to again deter- 
mine the dipole magnitude of the flight prototype thruster, as well as the actual flight 
thrusters and the overall S/C-SSU combination. The results of these determinations 
are also presented in table HI. It will be seen that the prototype value obtained by this 
method is about 22 percent higher than that obtained using the torsional pendulum setup. 

It is felt that the field measurement technique produces the more accurate value, and 
that the difference can be mainly attributed to the nonlinearity of the magnetic field 
lines for points outside the coil plane. In any case, the values are sufficiently close 
that assurance was gained as to the adequacy of the dipole compensation previously de- 
scribed. 

Measurements were made on the spacecraft- SSU combination in both the unpowered 

and powered configurations, to determine the contribution of any uncompensated current 

loops to the overall dipole. As can be seen from table m, these values are quite small 

compared to those which arise from permanent magnetism. It can also be seen that 

compensation by the two thruster magnets in the roll axis is adequate, although not as 

good as originally expected. This also indicates that the field measurement technique 

produced a more accurate value for the thruster dipole. 

The values received from the Agena contractor for the Agena vehicle and solar ar- 
4 4 2 

ray dipoles were 6. 0x10 and 1. 5x10 ampere-cm , respectively. No polarities were 

provided with these numbers and so they were assumed to add in worst-case fashion. 

Finally, an estimate based on solar array circuit diagrams indicated a solar array di- 

3 2 

pole component of 4x10 ampere-cm along the pitch axis. When these numbers are 
added to those obtained for the S/C-SSU combination, the following (worst case) totals 
were obtained: 

38 



(1) Roll: -9. 8xl0 4 A-cm 2 

(2) Pitch: -8. 2xl0 4 A-cm 2 

(3) Yaw: -9. 8xl0 4 A-cm 2 

In summary, it can be stated that, while the methods just described did not provide 
the magnetic dipole values with great accuracy, the estimates obtained yielded assurance 
that magnetic disturbances would be relatively small and would not compromise the suc- 
cess of the mission. 
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APPENDIX C 


SERT II EQUATIONS OF MOTION 

It will be convenient for linearization to relate the reference axis system and the 
body (or principal) axis system in terms of Euler angles. The three sequential rota- 
tions required to find the body axes are shown in figure 4. The rotation about Z r is a 
yaw angle i p; the rotation about y' is a pitch angle 9; and the rotation about x" is a 
roll angle ep. The reference and body systems are related by 

’x] op ce sp ce 

y = Ci p S9 Sep - Sip Cep Ci p Cep + S 9 Sip Sep 

z Ci p S9 Cep + Sip Sep Sip S9 Cep - Ci p Sep 

The inertial body angular rates can be expressed in terms of Euler angles, Euler 
angle rates, and the rotation rate of the reference axis system (orbit rate). 

u) = ep - ip sin 9 - co sin i p cos 9 (C2a) 

X LJ 

aiy = 9 cos ep + \p cos 9 sin ep - w Q ( cos ip cos ep + sin i p sin 9 sin ep) (C2b) 

• • 

u> z = -9 sin ep + ip cos 9 cos ep - o> 0 (-cos ip sin ep + sin ip sin 9 cos ep) (C2c) 

The satellite dynamics are expressed by Euler's equations of motion: 

W>x + u y u z^zz - V = T x < c3a ) 

+ Vz^x “ Kt) ~ Ty (C3b) 

^zz^z + ^x^y^yy ” ^xx^ - ^ z (C3c) 

T_ 

where Ijj is a principal moment of inertia and Tj is an external torque about the j in 
axis. 
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Gravity-Gradient Torques (Ref. 9) 


Assuming a body in earth orbit is attracted by an inverse-square gravitational force 
field, the force on a differential mass dm due to the external gravitational field is 


dF = - -U- r dm = 
r 3 




(C4) 


where 

H product of gravitational constant and earth mass 
r radius vector from earth center to dm 

r c radius vector from earth center to satellite center of mass 

p' vector locating dm with respect to the center of mass 

and r, r c , p’ are shown in figure 14. 

The differential torque about the center of mass is 

dT__ = p’ X dF 
gg F 

and the total gravity-gradient torque is 



Now since 

p’ x (r c + p') = P’ x r c 


(C 5) 


(C6) 


(C7) 
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equation (C6) becomes 


T = - JL 
gg 3 
r c 


/ 


i- 


3r • p’\ _ _ 


| (p f x r c )dm 


(C8) 


From figure 14, we note that in the i, j, k system 


p* = x i + yj + zk 


r c = r c (ai + bj + ck) 


Now, substituting these expressions into equation (C8) gives 


T = - JL 
gg 2 

r c 


/t 


1 _ 3(ax + by + cz) 


|(cy - bz)i + (ax - cx)j + (bx - ay)kjdm 


(C9) 


Now choosing i, j, k to be the principal axes of the body, the following relations hold: 

JL xdm =X y dm = zdm -° 

j£ xz dm =_^"yz dm = 0 


Z <y 2 - z2,dm = T zz - Vy 

f m (z 2 - x 2 )dm - ^ - I zz 

J t [ ( * 2 - y2)dm = Vy ' ! xx 
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I 




Thus equation (C9) becomes 



(CIO) 


Since a, b, c are the direction cosines of r c in the principal axis system, using 
equation (Cl) and remembering that r c is along the " Z R axis we find 


a = sin 9 


b = -cos 9 sin cp 
c = -cos 9 cos cp 

and 

2 

be = cos 9 sin cp cos < p 
ac = -sin 9 cos 9 cos cp 
ab = -sin 9 cos 9 sin cp 

Euler angles as 

(Cll) 

For an eccentric orbit, r c can be expressed as 

r - a(l - e 2 ) 
c 1 + e cos £ 

where 

e orbit eccentricity 


Thus the gravity-gradient torque may be expressed in terms of 

(l zz " l y^ cos2e sin V cos 

^gg = -(*xx ~ r zz> sin 0 cos 0 cos * 

r't 

-(I yy - ^ xx ) s i n 9 cos 9 sin cp 
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a semimajor axis 
£ true anomaly 

O 

For the SERT II orbit, e is small, so the quantity ju/rj; can be approximated as 


JL 


r 


3 

c 



(1 + 3e cos 



where 



and the time of perigee passage t is equal to 0. 

Control Moment Gyros 

As shown in figure 5, the CMG reference axes and the body (principal) axes can be 
related by 


’ORAj" 


'0 

0 

0 

X 


SRA i 

= 

0 

-cos cp gi 

-sin <p gi 

y 

(Cl 2) 

IRAj 


0 

sin <p gi 

-cos <p g . 

z 



where 


<p g l = +28° (CMG 1) 
cp g2 = -28° (CMG 2) 


Since the CMG’s are single-degree-of-freedom gyroscopes, the instantaneous true 
spin and input axes are related to the CMG reference axes by 
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0 



0 

COS Oj 

-sin cn 


sin 


COS Cn 


ORAj 


SRA, 


IRAj 


(C13) 


"f V*! 

where cn is the CMG gimbal angle of the i cn CMG and is positive for a positive ro- 
tation about the output axis. 

The CMG equations of motion may be expressed as 


Vi + Vi 


Vi = w IAi h i + w Di h i + T 


TGi 


(Cl 4) 


where 

gimbal inertia of V CMG 
c^ damping coefficient of i Ln CMG 

spring constant of V CMG 
h^ angular momentum of V CMG 

■Hi 

ojjAi true input axis rate of i CMG 
drift rate of V CMG 

r TGi programmed torque generator torque of V CMG 
Equation (C14) may be approximated by 


Vi + Vi = w IAi h i + T TGi < C15 > 

In equation (C15) we have neglected the CMG inertia and the drift rate. The inertia 
adds a negligible phase lag to the CMG response, and the drift rate in near zero grav- 
ity may be neglected. 

Now using equations (C12) and (C13) we can solve (C15) for d^: 


K' a ’‘ T TGi 


°i = ~ + 


c i c i 


h. 


+ — (wy( c °s cp gi sin cr. + sin <p gi cos cn) + w z (sin <p g . sin a. - cos <p gi cos 
c . L 

(C16) 
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Considering the CMG's as sources of external torque on the satellite, the torque 
on the body from the CMG’s may be expressed in the body axis system as 


T CMGi = " h i - w X 


(C17) 


Using equations (C12) and (C13), we find the CMG angular momentum in the body 
system 


h 



0 


l-cos <p gi cos a. + sin <p gi sin o\ 
-cos cn sin cp . - sin cos cp g . 


so 


■h = Vi 


-cos <pg i sin o\ - sin <p gi cos cr 
-sin sin a i + cos cp gi cos cn 


and 


(CIS) 


-a) xh 



(-C<p g . Co i + S<p gi SCT i ) - W y (-CCT. S<p gi - 
< jo (-Cc t. Sep . - Sct. Cep ■) 

X V 1 ^gl 1 ^gl ; 

-co x (- c <pgi Corj + S^ gi Sa i ) 


Summing equations (C18) and (C19) gives the CMG torques: 



(C19) 


T CMGi 



w z (-Cc p gi Ca t + Sc p gi Sct.) - co y (-Ccr. Scp gi - Sct. Ccp gi ) 
CT i (-Ccp gi Sct. - Sc p gi Ccr i ) + w x (-Cct. Scp gi - SOj C<p gi ) 
V _S< ?gi Sct. + Ccp gi CCTj) - w x (-Ccp gi Ccr i + Scp gi Sct.) 


(C20) 


CMG nonlinearities such as stiction, torquer fade, and drift rate are small (as 
shown in table n) and may be neglected in this analysis. These nonlinearities are in- 
cluded in the simulations. 
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Linearized Equations of Motion 


The linearized equations of motion are found by assuming that body and CMG 
angles and rates are small. Equations (C2) may be linearized as follows: 


w x = V - w o^ 

(C 21a) 

"y = 6 - u o 

(C21b) 

+ u Q <p 

(C21c) 

where 


w o = “o^ 1 + 2e cos w o^ 


Linearizing the derivatives of (C21) gives 


~ w o^ 

(C22a) 

co = 0 - co 
y o 

(C22b) 

W z = $ + U* 0 cp 

(C22c) 


where 


= -2e cu* 2 sin a;*t 
o o o 

Equations (Cll) may be linearized for small eccentricity: 


T ggx = - 3 ("o) 2 ( I yy - 

(C 23a) 

T ggy ' - 3 K) 2 «x X - 

(C 23b) 

r ggz = ' 3( "o> ^yy " 

(C23c) 
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If it is assumed that hj = h 2 = h, that Cj = C 2 = C, and that or. is small, then 
using equations (C20), the total linear CMG torques are 


T CMGx = ~ h j 2 *^ cos 28 ° + 2a) o^ cos 28 ° + w o cos 28 °( cr i + ct 2^ 

| (C24a) 

T CMGy = " h (°i “ CT 2^ sin 28 ° 

(C24b) 

T CMGz = + °’2^ cos 28 ° + 2h< ? cos 28 ° " cos 28° 

(C24c) 

Using equations (C16), we find 


= ~ ~ [^ 0(^1 ^ 2 ^®® 28 + 20 1 cos 28° + 2o) ^cp cos 28^*J 

(C25a) 

^1 “ °2 = ~ [ 2 ^ s * n 28 ° ” w o^ a l " ct 2^ cos 28 ° " 4ew o s * n 28 ° cos 

t] (C25b) 


where the torquer torque is 

T TGi = -K Sin ^gi 

Substituting equations (C21), (C22), (C23), and (C24) into equations (C3) gives 

W# ‘ ^ + ~ "o^ + V^zz ' V + 3(a, o )2(I yy " W* 

+ h^2i// cos 28° + 2co*(p cos 28° + + a 2 )cos 28°J = T x (C26a) 

i yy (e - w Q ) + W - « 0 U/)(* + ^ 0 <p)(i xx - i zz ) + 3(<) 2 (i xx - i zz )e 

+ Mdj - a 2 )sin 28° = T y (C26b) 

r zz^ + + (v " <o o'l / )0 - «o)(Iyy - r xx) + 2hu)*i// cos 28° 

- h(a 1 + a 2 )cos 28° - 2h<p cos 28° = T z (C26c) 
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The second terms of equations (C26a) and (C26c) give rise to the Euler cross-coupling 
terms. When linearized, the expressions are similar to the gravity-gradient torque 
terms. Considering the cross-coupling torques as external, the linearized forms are 


T x - -(“PVlyy - W 
T y = 0 

T x - -K^Oyy - W 

Equations (C26) maybe linearized further and together with equations (C25) completely 
describe the vehicle motion: 

l xx i + <p - I zz ) + 2hw* COB 28 °] + + [o,*(I yy - I xx - I zz ) + 2h cos 28°] 

+ hcu*(a^ + o^cos 28° = T x 

I yy 8 + 3(u)p 2 (I J[x - I zz )9 + 2e(ai*) 2 Iyy sin o*t + h(Oj - o 2 )sin 28° = T y 

hz* * - W + 2h “o cos 26 °] + i ["xx + l zz - V"o - 2h cos 28 °] 

- h(a^ + o^cos 28° = T z 

+ °2 = ~ “ [ w o (cr l + a 2^ cos 2 ®° + ^ cos 2 ®° + cos 


(C27a) 

(C27b) 

(C27c) 


CT 1 ” °2 = ~ ^ s ^ n 28° “ ^ a i - o'2^ cos 2 ®° “ 4eco* s * n 2 ®° cos <t] 

Neglecting initial conditions, the Laplace transformations of equations (C27) can be 
expressed as 
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(C28) 


m 

“ 




V* * 3 K> 2 «*x - W 

hs sin 28° 

9 (s) 


O)* 

T (s) - 2e(w*) 2 I 0 

y ' ' o> yy 2 , *,2 

s + {w* 0 ) 

-2hs sin 28° 

cs + hco* cos 28° 

cr'(s) 


-4heco* sin 28° s 





s 2 + (w *) 2 _ 


and 


W 2 + 4 “*ohyy - W + 2h < cos 28 ° 

sfcttyy • X xx - W + 2h cos 28 °] 

ha>* cos 28° 


<p{ s) 


T X (S) 

+ ! 2Z - V - 2h COS 28°] 

I zz s 2 + 2h U * 0 cos 26° + K) 2 (Iyy - ^ 

-hs cos 28° 

► 

^(s) 


T z (s) 

2hu>* cos 28° 
l o 

2hs cos 28° 

cs + hco* cos 28° 


_a M (s) 


0 


where 

°" = °1 - °2 
a ” = °1 + (J 2 


(C29) 
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APPENDIX D 





SATELLITE SIMULATION 

The satellite dynamics were simulated on both analog and digital computers. 


Analog 

The analog simulation involved rotational dynamics only. Two EAI 231 R-V com- 
puters were used for the simulation, which included all torques except magnetic torques. 
Solar pressure and eccentricity effects could be included without knowledge of the actual 
spacecraft position. 

The analog program was fast and easily adaptable. The studies performed using 
this program included general stability studies, satellite response to ion thruster oper- 
ation, and sensitivity to CMG parameters. 


Digital 

The digital simulation was performed on an IBM 7094-7044 Direct Couple system. 
The integration scheme was a fourth-order Runge-Kutta method including variable step- 
size and error control. The simulation included orbit dynamics and was used to study 
satellite response to magnetic torques. The program was also used to verify analog 
program results. 

Sketch (b) shows a general block diagram applicable to both analog and digital sim- 
ulations. The digital simulation includes all blocks shown. The analog simulation does 
not include satellite position simulation or magnetic torques. 
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APPENDIX E 


SYMBOLS 

A area, (ft^) 

a semimajor axis, m (ft) 

a, b, c direction cosines of r c 

C cosine 

c. damping constant of i th CMG, (N)(m)(sec) 

e eccentricity 

F force, N (lbf) 

H earth magnetic field intensity vector, Oe 

H e horizontal component of H, Oe 

h’ length of Agena/spacecraft/ spacecraft support unit, m (ft) 

In angular momentum of i Ln CMG, (N)(m)(sec) 

2 

I moment of inertia, kg-m 

IA- true input axis of i** 1 CMG 

IRAj input reference axis of i^ 1 CMG 

o 

Ijj satellite moments of inertia: i = x, y, z; j = x, y, z, kg-m 

i inclination, rad (deg) 

th 9 

J. gimbal inertia of i CMG, kg-m 

K- spring constant of i^ CMG, N-m/rad 

K w torsional spring constant, N-m/rad 

l component of thruster location along x axis, m 

l component of thruster location along z axis, m 

m mass, kg 

n unit vector normal to area 

OA i true output axis of i^ 1 CMG 

ORA^ output reference axis of i^ 1 CMG 

P g solar pressure, N/m^ (lbf/ft^) 

R reflectivity 



R e radius of earth, m 

r' Agena radius, m (ft) 

r radius vector from earth center to satellite differential mass dm, m (ft) 

r-g moment arm from satellite center of mass to point of force application, m (ft) 

r radius vector from earth center to satellite center of mass, m (ft) 

S sine 

SA i true spin axis of i CMG 
SRA i true reference spin axis of i in CMG 
s unit vector from satellite to sun 

T torque, N-m (ft-lbf) 

t time, sec 

a thruster pitch misalinement, rad (deg) 

a ' angle between -j axis and dA in i-j plane, rad (deg) 

/3 thruster roll misalinement, rad (deg) 

/3' instantaneous apparent sun declination referenced to satellite, rad (deg) 

y phase angle referencing the reference and Xg, Yg, Zg axis systems in time, 

rad (deg) 

5 apparent sun declination in the Xg, Yg, Zg axis system, rad (deg) 

£ true anomaly, rad (deg) 

77 angle between earth magnetic dipole axis and radius vector to satellite, rad 

(deg) 

0 pitch angle, rad (deg) 

3 2 

/i product of universal gravitation constant and earth mass, m /sec 

_ 0 

p satellite magnetic moment vector, A- cm 

— 2 

prj, thruster magnetic moment vector, A-cm 

p' vector location of differential mass dm with respect to satellite center of 

mass, m (ft) 

gimbal angle of i th CMG, rad (deg) 

it. 

r TGi programmed torquer torque of i tn CMG, N-m 
(p roll torque 
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^gi 

1L 

mounting angle of i n CMG, rad (deg) 



* 

yaw angle, rad (deg) 



n 

inertial longitude of orbit ascending node, rad (deg) 



n E 

inertial longitude of Greenwich, rad (deg) 



O? 

inertial body rate vector, rad/sec (deg/sec) 



U) 

0 

instantaneous orbit rate, rad/sec (deg/sec) 



CO* 

o 

average orbit rate, rad/sec (deg/sec) 



w Di 

drift rate of i** 1 CMG, rad/sec (deg/sec) 



w IAi 

■f-Vi 

input axis rate of i in CMG, rad/sec (deg/sec) 



Subscripts: 




A 

aft 



ax 

axial magnetic component 



B 

body (principal) axis system 



C 

coil 



CMG 

control moment gyro 



cc 

cross coupling 



F 

fore 



G 

geometric axis system 



GX, GY, GZ 

geometric axis system components 



gg 

gravity gradient 



I 

inertial axis system 



M 

magnetic axis system 



max 

maximum 



norm 

normal magnetic component 



O 

orbit axis system 



R 

reference axis system 



S 

sun axis system 



ss 

steady-state frequency response 



step 

steady-state step response 



T 

thruster axis system 



x, y, z 

body (principal)-axis-system components 






55 







REFERENCES 


1. Scott, E. D.: Control Moment Gyro Gravity Stabilization. Paper 63-324, AIAA, 

Aug. 1963. 

2. Fischell, Robert E. : Passive Gravity-Gradient Stabilization for Earth Satellites. 

Torques and Attitude Sensing in Earth Satellites. S. Fred Singer, ed. , Academic 
Press, 1964, pp. 13-30. 

3. DeLisle, J. E.; Ogletree, E. G. ; and Hildebrant, B. M. : The Application of Gyro- 

stabilizers to Orbiting Vehicles. Torques and Attitude Sensing in Earth Satellites. 

S. Fred Singer, ed. , Academic Press, 1964, pp. 31-72. 

4. Roberson, Robert E. : Generalized Gravity-Gradient Torques. Torques and Attitude 

Sensing in Earth Satellites. S. Fred Singer, ed. , Academic Press, 1964, pp. 

73-82. 

5. Cake, James E. : A Fortran Code for Computing the Principle Mass Moments of 

Inertia of Composite Bodies. NASA TM X-1754, 1969. 

6. Kendrick, J. B. : TRW Space Data. TRW Systems Group, 1967, p. 76. 

7. McElvain, R. J. : Effects of Solar Radiation Pressure upon Satellite Attitude Control. 

Paper 1918-61, ARS, Aug. 1961. 

8. Kraft, J. D. : Guidance and Control. Space Flight Handbooks. Vol. I. Orbital 

Flight Handbook. Part 3 - Requirements. NASA SP-33, Pt. 3, 1963, pp. Xn- 
40, 41. 

9. Greenwood, Donald T. : Principles of Dynamics. Prentice-Hall, Inc., 1965. 



TABLE I. - SERT II SPACECRAFT AND ORBIT PARAMETERS 


o 

Inertia, kg-m : 

Roll, 

Pitch, Iy y 

Yaw > X zz 

Spacecraft mass, kg 

Semimajor axis, m 

Eccentricity 

Inclination, deg 

Ion thruster thrust, N 

Fore thruster location from center of mass, m: 

Along z axis 

Along x axis 

Aft thruster location from center of mass, m: 
Along z axis 

Along x axis 

2 2 

Solar array area, m (ft ) 

o 

Cross-sectional area of Agent (SSU and S/C), m 

Reflectivity of solar array 

Reflectivity of Agena 


8.5832X10 3 

10.9451X10 3 

2. 7469X10 3 

1.4349X10 3 

7. 378xl0 3 (expected) 

0. 005 (expected) 0. 0003 (measured) 
. 99. 1 (expected) 99. 135 (measured) 
2. 758X10" 2 

+2.847 

+0.536 

+2.847 

-0. 536 

17.65 (190) 

(ft 2 ) 9. 01 (97) 

0 . 2 

0.8 


TABLE H. - CONTROL-MOMENT-GYRO PARAMETERS 


Angular momentum, h, (N)(m)(sec) 4. 5 

Damping constant, c, (N)(m)(sec) 2. 25 to 4. 5 

Spring constant 3, 0. 0 

_3 

Torquer torque, N-m 2. 11x10 

Torque fade a , N-m/deg 8. 7x10 6 

Stiction, N-m 1. 63x10 ® 

Drift, N-m 3. 25xl0“ 5 

Mounting angles (spin reference axes location from negative pitch 

axis in yaw- pitch plane), <p , deg ±28 

Gimbal stops, deg from spin reference axes ±30 


a Torquer fade acts as a very light spring in the CMG. Since spring constant is essen- 
tially zero, CMG behaves as rate-integrating gyro with a large angular momentum. 



TABLE III. 


MAGNETIC DIPOLE MOMENTS 


(a) Thrusters 



Torsional 

pendulum 



Field meas- 
urements 

Magnetic dipole moment, 
A-cm 2 

Flight prototype 
Flight thruster 1 
Flight thruster 2 

22.7X10 4 

27. 7xl0 4 
27. 5 
26.6 


(b) Spacecraft and spacecraft 
support unit 



Permanent 

dipole 

Current loop 
dipole 

Magnetic dipole moment, 
A-cm 2 

Roll axis 
Pitch axis 
Yaw axis 

-2. 27xl0 4 
-1. 11 
-2. 33 

0. 10X10 4 
. 35 
. 25 





'-Solar array 



Figure 1. - General arrangement of SERT II spacecraft/Agena D. 










Figure 5. - Control- moment-gyro axes in body axis system. Gimbal stops are located *30° from spin reference 
axis. Gimbal angle of CMG is a rotation about output axis and is referenced to spin reference axis. 



CMG 2 gfm- CMG 1 gimbal angle, Pitch angle, 0, deg Roll angle 

bal angle, Oj, deg <p, deg 
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